Non-equilibrium spin polarization effect in spin-orbit coupling system and contacting 

leads 



o 
o 

> 

o 

00 



CD 



I 

O 
O 



> 
O 

in 
o 



I 

O 

o 



Yongjin Jiang^'El 

^Physics department of ZheJiang Normal University, Jinhua, Zhejiang Province, 321004, P-R-China 

We study theoretically the current-induced spin polarization effect in a two-terminal mesoscopic 
structure which is composed of a semiconductor two-dimensional electron gas (2DEG) bar with 
Rashba spin-orbit (SO) interaction and two attached ideal leads. The nonequilibrium spin density 
is calculated by solving the scattering wave functions explicitly within the ballistic transport regime. 
We found that for a Rashba SO system the electrical current can induce spin polarization in the SO 
system as well as in the ideal leads. The induced polarization in the 2DEG shows some qualitative 
features of the intrinsic spin Hall effect. On the other hand, the nonequilibrium spin density in the 
ideal leads, after being averaged in the transversal direction, is independent of the distance measured 
from the lead/SO system interface, except in the vicinity of the interface. Such a lead polarization 
effect can even be enhanced by the presence of weak impurity scattering in the SO system and may 
be detectable in real experiments. 

PACS numbers: 72.25.-b, 75.47.-m 



I. INTRODUCTION 



Recently there has been hot research interest in spin- 
polarized transport phenomena in semiconductor sys- 
tems with spin-orbit ( SO ) coupling because of their 
potential applications in the field of semiconductor 
spintronics^i^. One of the principal challenges in semi- 
conductor spintronics is the efficient injection and ef- 
fective manipulation of non-equilibrium spin densities 
and/or spin currents in nonmagnetic semiconductors by 
electric means, and within this context, the recently pre- 
dicted phenomenon of intrinsic spin Hall effect would 
be much attractive. This phenomenon was conceived 
to survive in some semiconductor systems with intrinsic 
spin-orbit interactions and it consists of a spin current 
contribution in the direction perpendicular to a driving 
charge current circulating through a sample. This phe- 
nomenon was firstly proposed to survive in p-doped bulk 
semiconductors with spin-orbit split band structures by 
Murakami et ali^ and in two-dimensional electronic sys- 
tems with Rashba spin-orbit coupling by Sinova et al.fS 
In the last two years lots of theoretical work have been 
devoted to the study of this phenomenon and intensive 
debates arose on some fundamental issues concerning this 
effect^i^^^^iii^ii^^i^^il^ia^ Now it was well es- 
tablished that this effect should be robust against im- 
purity scattering in a mesoscopic semiconductor system 
with intrinsic spin-orbit coupling, but it cannot survive 
in a diffusive two-dimensional electron gas with k-linear 
intrinsic spin-orbit coupling due to vertex correction from 
impurity scatteringi2iiSiiiii2£. It also should be noted 
that, prior to the prediction of this effect, a similar phe- 
nomenon called extrinsic spin Hall effect had also been 
predicted by Dyakanov-Perel^i and Hirsc h^^'^'^ . Unlike 
the intrinsic spin Hall effect, the extrinsic spin Hall ef- 
fect arises from spin-orbit dependent impurity scattering 
and hence it is a much weak effect compared with the 
intrinsic spin Hall effect, and especially, it will disappear 



completely in the absence of spin-orbit dependent impu- 
rity scattering. 

A physical quantity which is intuitively a direct con- 
sequence of a spin Hall effect in a semiconductor strip is 
the nonequilibrium spin accumulation near the transver- 
sal boundaries of the strip when a longitudinal charge 
current circulates through it. The recent experimental 
observations in n-doped bulk GaAs^l and in p-doped 
two-dimensional GaAs^ have revealed the possibility of 
generating non-equilibrium edge spin accumulation in 
a semiconductor strip by an extrinsio^ or intrinsic^— 
spin Hall effect, though from the theoretical points of 
view some significant controversies may still exist for 
the detailed theoretical interpretations of some experi- 
mental dataii^ii^i2ii2^i^. It also should be noticed that, 
in addition to the spin Hall effect, there may be some 
other physical reasons that may lead to the generation 
of electric-current-induced spin polarization in a spin- 
orbit coupled system. For example, the experimental re- 
sults reported in Re:2§ had demonstrated another kind of 
electric-current-induced spin polarization effect in a spin- 
orbit coupled system, which can be interpreted as the in- 
verse of the photogalvanic effect^^ . In the present paper 
we are interested in the question that, besides the trans- 
verse boundary spin accumulations, what kinds of other 
physical consequences can be predicted for a semicon- 
ductor strip with intrinsic SO coupling. To be specific, 
we will consider a strip of a mesoscopic two-dimensional 
electron gas with Rashba spin-orbit coupling connected 
to two ideal leads. For such a system, we show that when 
a longitudinal charge current circulates through the strip, 
besides the transverse boundary spin accumulations in 
the strip, non-equilibrium spin polarizations can also be 
induced in the leads, and the non-equilibrium spin polar- 
izations in the leads will also depend sensitively on the 
spin-orbit coupling parameters in the strip. Furthermore 
we will show that in the absence of spin-orbit coupling in 
the leads, the non-equilibrium spin polarizations in the 
leads will be independent of the distance away from the 
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border between the lead and the strip, except in the vicin- 
ity of the border where local states play an important 
role. The signs and the amplitudes of the non-equilibrium 
spin polarization in the leads is adjustable by switching 
the charge current direction and by tuning the Rashba 
SO coupling strength in the strip. Furthermore, in the 
presence of weak disorder inside the SO strip, the spin po- 
larization effect in the leads can be even enhanced while 
spin polarization in the strip is decreased, so such an 
effect might be more easily detectable than the spin po- 
larization inside the strip. The non-equilibrium lead po- 
larization effect provides a new possibility of all-electric 
control on spin degree of freedom in semiconductor and 
thus have practical application potential. 

The paper will be organized as following: In section 
II the model Hamiltonian and the theoretical formalism 
used in the calculations will be introduced. In section III 
some symmetry properties and numerical results about 
the non-equilibrium spin polarizations both in the strip 
and in the leads will be presented and discussed. At last, 
we summarize main points of the paper in section IV. 
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FIG. 1: schematic geometry of multi-terminal setup. 



II. METHOD OF SOLVING SCATTERING 
WAVE FUNCTION AND MULTI-TERMINAL 
SCATTERING MATRIX 

The system that we will consider in the present paper 
is depicted in Fig. [2 which consists of a strip of a meso- 
scopic Rashba two-dimensional electron gas connected to 
two half-infinite ideal leads. Each lead is connected to an 
electron reservoir at infinity which has a fixed chemical 
potential. In the tight binding representation, the Hamil- 
tonian for the total system reads: 

p <ij>a Ri.a 
<Ri,Rj>a 



Ri 



-t Y (CI^Cr,^, + h.c.) 

P,^,Rr^ 



(1) 



t 2 

Here t = 2m* a'^ hopping parameter between 

two nearest-neighbor sites where m* is effective mass of 



electrons and a the lattice constant in the 2DEG bar. 

= (C'fl.^i, Cfl.^|) is the spinor vector for the site Ri 
( Ri = {x,y)) in the Rashba SO system, a^^a^ are the 
standard pauli matrices. Cp^a is the annihilation oper- 
ator for the site pj with spin a in lead p. p„ and i?„ 
stand for nearest-neighbor pair on the two sides across 
the interface of the SO system and the lead p. wr^ is the 
on-site energy in the SO system, which can readily incor- 
porate random impurity potential. In a pure system we 
always set wr^ to zero. Ir is Rashba coupling coefficient. 

Our calculations will follow the same spirit of the usual 
Landauer-Biittiker's approach. But unlike the frequently 
used Green's function formalis m^^i*^" , in the formalism 
used in the present paper we will solve the scattering 
wave functions in the whole system explicitly. To this 
end, we firstly consider the scattering wave of an elec- 
tron incident from a lead. The real space wave function 
of an incident electron with spin a will be denoted as 
e~^'^^^i'Xm(y{yp)i where XmaiVv) denotes the m'th trans- 
verse mode with spin index a in the lead p and fc*^ the 
longitudinal wave vector. We adopt the local coordinate 
scheme for all leads. In the local coordinate scheme, 
the longitudinal coordinate Xq in lead q will take the 
integer numbers l,2,...,oo away from the 2DEG inter- 
face and the transverse coordinate Uq take the value of 
~Nq/2, ...,Nq/2. The longitudinal wave vector k^^ satis- 
fies the relation — 2icos(fcPj) + ef^ = E, where ef^ is the 
eigen-energy of the m'th transverse mode in lead p and 
E the energy of the incident electron. Including both the 
incident and reflected waves, the total wave function in 
the lead q has the the following general form: 



) = SpqS, 



■iVv 



,y1 



\vd (2) 



where (j^^^^t stands for the scattering amplitude from the 
(ma) mode in lead p to the (ncr') mode in lead q. The 
scattering amplitudes will be determined by solv- 

ing the Schrodinger equation for the entire system, which 
has now a lattice form and hence there is a separate equa- 
tion for each lattice site and spin index. Since Eq. ([2]) is a 
linear combination of all out-going modes with the same 
energy E, the Schrodinger equation is satisfied automat- 
ically in lead g, except for the lattice sites in the first 
row ( i.e., = 1 ) of the lead that are connected di- 
rectly to the 2DEG bar. The wave function in the first 
row of a lead, which are determined by the scattering 
amplitudes (j}^^^,, must be solved with the wave func- 
tion in the 2DEG bar simultaneously due to the coupling 
between the lead and the 2DEG bar. To simplify the no- 
tations, we define the wave function in the 2DEG bar as 
a column vector ip whose dimension is 2N ( TV is the total 
number of lattice sites in the 2DEG bar ). The scattering 
amplitudes (j)^^^^, will be arranged as a column vector (j) 
whose dimension is 2M ( M is the total number of lattice 
sites in the first row of the leads ). From the lattice form 
of the Schrodinger equations for the 2DEG bar and the 
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first row of a lead, one can obtain the following equations 
reflecting the mutual influence between the two parts: 

A^} = b + B0, 

C(j} = d + DV^, (3) 

where A and C are two square matrices with a dimension 
of 27V X 2N and 2M x 2M, respectively; B and D are 
two rectangular matrices describing the coupling between 
the leads and the 2DEG bar, whose matrix elements will 
depend on the actual form of the geometry of the system. 
The vectors b and d describe the contributions from the 
incident waves. Some details of deduction as well as the 
elements of these matrices and vectors have been given 
elsewhere^. 

After obtaining all scattering amplitudes i'^^a' 
, we can calculate the charge current in each 
lead through the Landauer-Buttiker formula, Ip = 
i^^lh) E.,... {T^al K-TP-^^Vp), where V, = nj{~e) 
is the voltage applied in the lead q and /i^ is the chemical 
potential in the lead g, T^^, are the transmission prob- 
abilities defined by = Y.m,n\^lna'?^qn/vpm and 
Vpm — 2tsm{kP^) is the velocity for the m'th mode in 
the lead p. 

With the wave function ^/'^7"^(i?i) in the 2DEG strip 
at hand, the non-equilibrium spin density in the 2DEG 
strip can also be calculated readily by taking proper en- 
semble average following Landauer's spirit--. We as- 
sume that the reservoirs connecting the leads at in- 
finity will feed one-way moving particles to the leads 
according to their own chemical potential. Following 
the Landauer's spirit^, let's normalize the scattering 
wave function tj^P^'^ so that there is one particle for 
each incident wave, i.e., we normalize e'^'^-^^i'Xmaiyp) 
to e~*'^™^''x^^(?/p)/vT , where i ^ oo is the length 
of lead p. Meanwhile, the density of states in lead p is 
l^m = ' where Wp„ = 2tsin(fcPj) is the veloc- 

ity of the m'th transverse mode in lead p. Adding the 
contributions of all incident channels of lead p with cor- 
responding density of states ( DOS ), we can obtain the 
total non-equilibrium spin density. In the linear trans- 
port regime, the spin density can be calculated with the 
incident energy at the Fermi surface as: 

(4) 

where {Sa{Ri)) denotes the spin density at a lattice site 
Ri in the 2DEG strip and /Zp is the chemical potential of 
lead p. 



III. RESULTS AND DISCUSSIONS 

When the system is in the equilibrium state, there will 
be no net spin density since the Hamiltonian has time- 
reversal { T ) symmetry. However, as a charge current 



flows from lead 1 to lead 2, non-equilibrium spin den- 
sity may emerge in the whole system, both inside the 
Rashba bar and in the leads. In this section we will 
investigate the non-equilibrium spin density under the 
condition of a fixed longitudinal charge current density. 
In our calculations we will take the typical values of the 
electron effective mass m* = O.OAme and the lattice con- 
stant a = 3nm'^4,. The chemical potential difference be- 
tween the two leads will be set by fixing the longitu- 
dinal charge current density to the experimental value ( 
w 100/iA/1.5/iTO ) as reported in Re^ and the Fermi en- 
ergy of the 2DEG bar will be set to Ef — — 3.8t through- 
out the calculations. We will limit our discussions to the 
linear transport regime at zero temperature. 
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FIG. 2: (Color online) Typical pattern of current-induced non- 
equilibrium spin density (in units of ^) in a two terminal 
structure. The lattice size for the central SO bar is 40 x 40. 
The Rashba coupling strength is taken as tR—O.lt. 

To show clearly the characteristics of the current- 
induced spin polarization in a two-terminal structure, be- 
low we will study the spin density in the Rashba bar and 
in the leads separately. Firstly we study the spin den- 
sity in the Rashba bar. The typical spin density pattern 
obtained in a two terminal structure is shown in Figj^l 
From the figures one can see that the spatial distribution 
of the spin density inside the Rashba bar exhibits some 
apparent symmetry properties, which are summarized in 
Eq.® given below. Theoretically speaking, these sym- 
metry properties are the results of some symmetry op- 
erations implicit in the two-terminal problem. Let us 
explain this point in some more detail. From the sym- 
metry property of the Hamiltonian ([1]) and the rectangu- 
lar geometry of the system shown in Figll] one can see 
that in our problem we have two symmetry operations 
icTyPx and iaxPy, where Px and Py denote the spacial re- 
flection manipulation y —y and x — > —x respectively 
and iax and iuy the spin rotation manipulation with the 
angle tt around the Sx and Sy axes in spin space respec- 
tively. By considering these symmetry operations, from 
Eq.(4) one can obtain immediately the following symme- 
try relations^: 

< Sx.,z{x,y) >i = ~ < Sx.z{x,-y) >i 
< Sy{x,y) >i = < Sy{x,~y) >i 
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<Sx{x,y)>i = -<Sx{-x,y)>i 
< Sy^z{x,y) >i = < Sy^^{-x,y) >i, (5) 

where < • • • >/ stands for the spin density induced 
by a longitudinal charge current / flowing from lead 1 
to lead 2. The first two lines of Eq.® result from 
the symmetry operation iayP^, which has been known 
widely before^'^^. The second two, however, are re- 
sults from the symmetry operation iUxPy and the T- 
reversal operation together—. Firstly, due to the T- 
reversal invariance, the equilibrium spin density vanishes 
when all chemical potentials are equal in Eq.(I3]), i.e., 

< Sa{x,y) >eq^< Sa{x,y) >/ + < Sa{x,y) >_/= 0, 

where < Saix,y) >-/ denotes the spin density induced 
by a longitudinal charge current / from lead 2 to lead 1. 
Secondly, due to the geometry symmetry of the system 
under the manipulation iaxPy, we have < Sx{x,y) >i—< 
Sxix,y) >-i and < Sy^z{x,y) >/= - < Sy^z{x,y) >_/. 
Combining these two results, we get the last two lines in 
Eq.®. 

In Fig[3]we plot the longitudinally averaged spin den- 
sity inside the Rashba bar as a function of the transverse 
position y, where < S{y) >^denotes the longitudinal av- 
eraged value of < S{x,y) > along the x direction ( i.e., 
the direction of the charge current flow ). Due to the 
symmetry relations shown in Eq.(IS]), the non- vanishing 
components are < Sy,z{y) >^ ■ Moreover, under the spa- 
tial reflection manipulation y — > — y, < Sy{y) >^ is 
even and < Sz{y) >^ is odd. The fact that the out- 
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FIG. 3: (Color online) Non-equilibrium spin density averaged 
in longitudinal direction in SO system , < Sy,z >^, as a 
function of transverse position, y. The y and z components 
has diflerent reflection symmetry with respect to x axis. The 
sign near the boundaries can be flipped by tuning Rashba 
coupling strength tR. 

of-plane component < Sz{y) >^ has opposite signs near 
the two lateral edges is consistent with the phenomenol- 
ogy of SH E^^i^^ . But for a two-terminal lattice structure 
with general lattice sizes, our results show that the spin 
density does not always develop peak structures near the 
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FIG. 4: (Color online)The variation of spin density averaged 
in transverse direction of lead 2 versus Rashba coupling con- 
stant. 



two boundaries but oscillate across the transverse direc- 
tion. This is somehow different from the naive picture 
of spin accumulation near the boundaries due to a spin 
Hall current. The in- plane spin polarization is not re- 
lated to the phenomenology of spin Hall effect. It can 
be regarded as a general magnetoelectric effect due to 
spin-orbit coupling"^^. It should be noted that, from the 
theoretical points of view, the relationship between spin 
current the induced spin polarization is actually a much 
subtle issue and is currently still under intensive debates. 
In this paper, however, we will free our discussions from 
such controversial issues. 

It is evident from Fig. [3] that not only the magni- 
tude, but also the sign of the spin accumulation near 
the two transverse boundaries can be changed by tun- 
ing the Rashba coupling strength. This is also reported 
in Re^ for a continuum two-terminal model. Since the 
Rashba coupling strength can be tuned experimentally 
through gate voltage, such a flipping behavior for spin 
accumulation might provide an interesting technological 
possibility for the electric control of the spin degree of 
freedom. 

Now we present the most important theoretical predic- 
tion of this paper, i.e., the charge-current induced spin 
polarization effect in the contacted leads. Due to the first 
two lines in Eq. ^ , after taking average in the transversal 
direction ( which will be denoted by < Sa >^ as a func- 
tion of the longitudinal coordinate X2 ), only < Sy >^ 
will remain nonzero. In Fig.[4]< Sy >^ is plotted versus 
the longitudinal coordinate x of the whole system. The 
lattice size of the whole system is set to 120 x 40, while the 
central 40 x 40 lattice sites represents the Rashba SO bar 
(see FigH]). As indicated clearly in the figure, the border 
between the Rashba bar and the contacted leads are at 
X = ±20. The Rashba SO coupling strength is chosen 
to be tBjt = 0.05, 0.1, 0.15. As can be seen from Fig. gl 
< Sy >^ changes with the coordinate x both inside the 
Rashba bar and in the leads, but in the leads it changes 
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with the coordinate x only in a much narrow region close 
to the border between the leads and the Rashba bar and 
will reach to a fixed value further away from the border. 
This phenomenon is the lead spin polarization effect in 
this paper. Theoretically, this phenomenon can be un- 
derstood as following. For an incident wave (m, a) with 
a Fermi energy Ef in lead 1, the scattering wave func- 
tion in lead 2 is: Vi^'^(a;2, yz) = E„ &e*'^""=Xn(y2), 
where — 2tcos(fc„) + £„ = Ef. Then the spin density 

< S{x2, 1/2) > can be calculated as ( dropping off a nor- 
malization factor ), 

E — E ^"""(^2, 2;2)aa,/3V'r"(x2, ys), 

Thus, after taking average along the transversal direc- 
tion, wc have 

< S{x2) (l/iV2)E < S{x2,y2) > 
y2 

2^[02nae Cra,/3(</>2n/3 e ')J(6) 

m,<j 1™ n,a,[3 

where we have used the orthogonality relation for the 
transverse modes: I^yj X«(y2)x«'(y2) = 5«n'- Evidently, 
the summand [• • •] in Eq. ([6]) is independent of X2 as long 
as fc„ is real ( i.e., for those longitudinal propagating 
modes ) . If fc„ is imaginary, the corresponding mode will 
describe an exponentially localized state in lead 2 in the 
vicinity of the border between the lead and the Rashba 
bar. Such evanescent components^! can contribute to the 
local charge density and spin density only in the vicin- 
ity of the border. At some distance far away from the 
interface, the contribution of the evanescent modes to 

< S{x2) >"^will decay exponentially as X2 increase and 
hence < S{x2) >^ will be independent of X2, as illus- 
trated in Fig. m This result implies that the spin po- 
larization inside the Rashba bar can be inducted out by 
the driving charge current to the leads and manifest itself 
in an amplified way, i.e., the leads will become spin po- 
larized and thus an electric-controllable spin state in the 
leads can be realized. Such a lead spin polarization effect 
might be more easily detectable by magnetic or optical 
methods than the spin polarization inside the mesoscopic 
Rashba bar. 

The results presented above is obtained in the absence 
of impurity scattering. To simulate spin-independent dis- 
order scattering, we assume a uniformly distributed ran- 
dom potential w € [—Wd, Wd] in the Rashba bar with a 
disorder strength Wd ■ The non-equilibrium spin density 
can be obtained by taking average for a number of ran- 
dom realizations of local potentials. We averaged 1000 
random realizations in all calculations. In Fig. [5] we show 
how the curves of < S*;^ >^ versus x change with Wd- 
All curves are calculated with t^ = Q.lt. Noticeably, in- 
side the Rashba bar the height of < Sy >"^decreases with 
increasing Wd but the lead spin polarization will firstly 
increase in the weak disorder regime ( Wd < 0-5t ) and 



FIG. 5: (Color online)The variation of spin density averaged 
in transverse direction versus longitudinal coordinate at the 
presence of disorder scattering. 



then decrease in the stronger disorder regime ( note that 
our calculations are carried out under the condition of 
a fixed longitudinal charge current density). Thus, there 
exists a certain regime in which disorder can enhance the 
saturation value of the spin polarization in the leads. The 
slight asymmetric form of Wd — '2t line means that we 
need to take more random configurations of local poten- 
tial for large Wd value, which is reasonable. Of course, 
there are some other factors such as spin-orbit coupling 
or impurity scattering in the leads that may reduce or 
suppress the lead spin polarization effect. Further quan- 
titative study is needed in order to clarify the role of 
these factors. 



IV. CONCLUSION 

To conclude, in this paper we have presented a the- 
oretical study on the non-equilibrium lead spin polar- 
ization effect in a two-terminal mesoscopic Rashba bar 
under the condition of a fixed longitudinal charge cur- 
rent density. We have predicted that a finite amount of 
non-equilibrium spin polarizations can be induced in the 
leads by the spin-orbit coupling inside the mesoscopic 
Rashba bar when the longitudinal charge current cir- 
culate through it. Such a lead spin polarization effect 
can survive in the presence of weak disorder inside the 
Rashba bar and thus might be observable in real ex- 
periment. Such an effect might provide a new kind of 
electric-controllable spin state which is technically attrac- 
tive. But it should be stressed that, for real systems, due 
to the existence of disorder scattering or spin-orbit inter- 
actions or other spin decoherence effects in the leads, the 
lead spin polarization effect predicted in the present pa- 
per may be weakened. These factors need to be clarified 
by more detailed theoretical investigations in the future. 
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